Compact star constraints on the high-density EoS 
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Abstract. A new scheme for testing the nuclear matter 
(NM) equation of state (EoS) at high densities using con- 
straints from compact star (CS) phenomenology is applied 
to neutron stars with a core of deconfined quark matter 
(QM). An acceptable EoS shall not to be in conflict with 
the mass measurement of 2.1 ± 0.2 Mq (1 a level) for 
PSR J0751+1807 and the mass-radius relation deduced 
from the thermal emission of RX J1856-3754. Further con- 
straints for the state of matter in CS interiors come from 
temperature-age data for young, nearby objects. The CS 
cooling theory shall agree not only with these data, but 
also with the mass distribution inferred via population 
synthesis models as well as with LogN-LogS data. The 
scheme is applied to a set of hybrid EsoS with a phase 
transition to stiff, color superconducting QM which fulfills 
all above constraints and is constrained otherwise from 
NM saturation properties and flow data of heavy-ion colli- 
sions. We extrapolate our description to low temperatures 
and draw conclusions for the QCD phase diagram to be 
explored in heavy- ion collision experiments. 



1. Introduction 

Recently, new observational limits for the mass and 
the mass-radius relationship of CSs have been obtained 
which provide stringent constraints on the equation of 
state of strongly interacting matter at high densities, see 



Klahn et al. (2006) and references therein. In this latter 



work several modern nuclear EsoS have been tested re- 
garding their compatibility with phenomenology. It turned 
out that none of these nuclear EsoS meets all constraints 
whereas every constraint could have been fulfilled by some 
EsoS. As we will point out in this contribution, a phase 
transition to quark matter in the interior of CSs might 
resolve this problem. In the following we will apply an ex- 
emplary EoS for NM obtained from the ab-initio relativis- 
tic Dirac-Brueckner-Hartree-Fock (DBHF) approach us- 



ing the Bonn A potential (van Dalen et al. (2005) ). There 



it would require an essentially nonperturbative treatment 
of QCD at finite chemical potentials. For some promis- 
ing steps in the direction of a unified QM-NM description 
on the quark level, we refer to the nonrelativistic poten- 



tial model approach by Ropke et al. (1986) and the NJL 
model one by Lawley et al. (2006) Simulations of QCD on 
the Lattice meet serious problems in the low-temperature 
finite-density domain of the QCD phase diagram relevant 
for CS studies. However, there are modern effective ap- 
proaches to high-density QM which, albeit still simplified, 
focus on specific nonperturbative aspects of QCD. They 
differ from the traditional bag model approach and al- 
low for CS configurations with sufficiently large masses. 
For our QM description we em- 



Alford et al. (2006) 



see 

ploy a three-flavor chiral quark model of the NJL type 
with selfconsistent mean fields in the scalar meson (cou- 
pling Gs) and scalar diquark (coupling Go ~ Vd Gs) 
channels ( Blaschke et al. (2005)j ), generalized by includ- 
ing a vector meson mean field (coupling Gv = W Gs), 
see 



Klahn et al. (2006a) 



We show that the presence of a QM core in the in- 
terior of CSs does not contradict any of the discussed 
constraints. Moreover, CSs with a QM interior would be 
assigned to the fast coolers in the CS temperature-age di- 
agram. Another interesting outcome of our investigations 
is the prediction of a small latent heat for the deconfine- 
ment phase transition in both, symmetric and asymmetric 
NM. Such a behavior leads to hybrid stars that "masquer- 
ade" as neutron stars and has been discussed earlier by 



Alford et al. (2005)| for a different EoS. This finding is of 
relevance for future heavy-ion collision programs at FAIR 
Darmstadt. 

2. The flow constraint from HICs 

The behaviour of elliptic flow in heavy-ion collisions 
is related to the EoS of isospin symmetric matter. 
The upper and lower limits for the stiffness deduced 



is not yet an ab-initio approach to the high-density EoS 
formulated in quark and gluon degrees of freedom, since 



from such analyses (Danielewicz et al. (2002)) are in- 
dicated in Fig. 1 as a shaded region. The nuclear 
DBHF EoS is soft at moderate densities with a com- 
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Fig. 1. Constraint on the high-density behavior of the 
EoS from simulations of flow data from heavy-ion collision 



experiments (shaded area from Danielewicz et al. (2002) ) 
compared to the nuclear matter and hybrid EsoS discussed 
in the text. 



pressibility K = 230 M eV ( |van Dalen et al. (2004)| 



Gross-Boelting et al. (1999) I, but tends to violate the flow 
constraint for densities above 2-3 times nuclear saturation. 
As a possible solution to this problem we adopt a phase 
transition to QM with an EoS fixed to sketch the upper 
boundary of the flow constraint. In order to obtain an EoS 
as stiff as possible we use a vector coupling of -qv — 0.50 
and a diquark coupling oi rjD — 1.03. Herewith the EoS is 
completely fixed. 

3. Constraints from astrophysics 

3.1. Maximum mass and mass-radius constraints 

These most severe constraints come in particular 
from the mass measurement for PSR J0751-I-1807 



(Nice et al. (2005) I giving a lower limit for the maximum 
mass « 1.9 Mq at la level, and from the thermal emis- 



sion of RX J1856-3754 (Triimper et al. (2004)) providing 



a lower limit in the mass-radius plane with minimal radii 
i? > 12 km. These constraints can only be fulfilled by a 
rather stiff EoS. The most stiff quark matter contribution 
to the EoS which still fulfills the flow constraint in sym- 
metric matter corresponds to r]v — 0.5 with a maximum 
mass for hybrid stars « 2.1 Mq, rather independent of 
the choice of r]o which fixes the critical mass for the on- 
set of deconfinement, see Figs. 2, 3. For a more detailed 
discussion, see Klahn et al. (2006)] Klahn et al. (2006a) 



3.2. Cooling constraints 

Direct Urea (DU) processes are flavor-changing pro- 
cesses with the prototype being n — > p + e~ -\- 
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Fig. 2. Stable CS conflgurations for neutron stars 
(DBHF) and hybrid stars, characterized by the param- 
eters r]D and r]v of the quark matter EoS. 
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Fig. 3. Mass-radius relations for CSs with possi- 
ble phase transition to deconfined quark matter, see 



Klahn et al. (2006a) 



(Gamow and Schoenberg (1941) ), providing the most 



effective cooling mechanism in the hadronic layer of 
compact stars. It acts if the proton fraction x exceeds 
the DU threshold xdUi ^ — np/{nn + rip) > xdu- The 



threshold is given by xdu — 0.11 (Lattimer et al. (1991) I 
and rises up to xdu — 0.14 upon inclusion of muons. 
Although the onset of the DU process entails a sensi- 
ble dependence of cooling curves on the star masses, 
hadronic cooling with realistic pairing gaps is not suffi- 
cient to explain young, nearby X-ray dim objects, like 
Vela, with typical CS masses, not exceeding 1.5 M© 
([Blaschke et al. (2004)[ [Grigorian et al. (2005)[). The 



point on the stability curve in Fig. 2 marks the DU 
threshold density for the DBHF EoS. Quark matter 
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Fig. 4. Cooling evolution for hybrid stars of different 
masses given in units of M©. Note that Vela is described 
with a typical CS mass not exceeding 1.45 Mq. From 
Popov et al. (2006a)[ 
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Fig. 5. Comparison of observational data for the 
LogN-LogS distribution with results from popula- 
tion synthesis using hybrid star cooling according to 
Popov et al. (2006a)[ 



DU processes provide enhanced cooling 
by the diquark pairing gaps 

Page et al. (2000) ) and their density dependence 
( [Grigorian et al. (2005)[ [Popov et al. (2006a)[ ). For a 



characterized 
( Blaschke et al. (2000)] 



Sedrakian (2007) 



recent review, see 

To verify this rather heuristic approach we apply ex- 
plicit calculations of the cooling of hybrid configurations 
which shall describe present data of the temperature- age 
distribution of CSs. The main processes in nuclear matter 
that we accounted for are the direct Urea, the medium 
modified Urea and the pair breaking and formation pro- 
cesses. Furthermore we accounted for the IS'o neutron 
and proton gaps and the suppression of the ZP2 neutron 
gap. For the calculation of the cooling of the quark core 
we incorporated the most efficient processes, namely the 
quark modified Urea process, the quark bremsstrahlung, 
the electron bremsstrahlung and the massive gluon-photon 
decay. In the 2-flavor superconducting phase one color of 
quarks remains unpaired. Here we assume a small residual 
pairing (Ax) of the hitherto unpaired quarks. For detailed 
discussions of cooling calculations and the required ingre- 
dients see [Blaschke et al. (2004)[ [Popov et al. (2006a)[ [D 
and references therein. The resulting temperature-age re- 
lations for the introduced hybrid EoS are shown in Fig. 4. 
The critical density for the transition from nuclear to 
quark matter has been set to a corresponding CS mass of 
Merit = 1.22 Mq. All cooling data points are covered and 
correspond to CS configurations with reasonable masses. 
In this picture slow coolers correspond to light, pure neu- 
tron stars (M < Merit), whereas fast coolers are rather 
massive CSs {M > Merit) with a QM core. 

Another constraint on the temperature-age relation is 
given by the maximum brightness of CSs, as discussed by 



Grigorian (2006) It is based on the fact that despite many 
observational efforts one has not observed very hot NSs 
(logT > 6.3 - 6.4 K) with ages of 10^ - 10^-^ years. Since 
it would be very easy to find them - if they exist in the 
galaxy - one has to conclude that at least their fraction is 
very small. Therefore a realistic model should not predict 
CSs with typical masses at temperatures noticeable higher 
than the observed ones. The region of avoidance is the 
hatched trapezoidal region in Fig. 4. 

The final CS cooling constraint in our scheme is 
given by the Log N-Log S distribution, where N is 
the number of sources with observed fluxes larger than 
S. This integral distribution grows towards lower fluxes 
and is inferred, e.g., from the ROSAT all-sky survey 



(Neuhauser and Triimper (1999)). The observed Log N- 



Log S distribution is compared with the ones calculated 
in the framework of a population synthesis approach in 
Fig. 5. A detailed discussion of merits and drawbacks can 



be found in Popov et al. (2006) 



Altogether, the hybrid star cooling behavior obtained 
for our EoS fits all of the sketched constraints under the 
assumption of the existence of a 2SC phase with X-gaps. 

4. Outlook: The QCD phase-diagram 

Within the previous sections we exemplified how to ap- 



ply the testing scheme introduced in Klahn et al. (2006) 



to the modeling of a reliable hybrid EoS with a NM-QM 
phase transition that fulfills a wide range of constraints 
from HICs and astrophysics. In a next step we extend 
the description to finite temperatures focusing on the be- 
haviour at the transition line. For this purpose we ap- 
ply a relativistic mean-field model with density-dependent 



masses and couplings (Typel (2005)) adapted such as to 
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Fig. 6. Phase diagram for isospin symmetry using the 
most favorable hybrid EoS of the present study. The NM- 
2SC phase transition is almost a crossover. The model 
DD-F4 is used as a finite-temperature extension of DBHF. 
For the parameter set (?7d = 0.75, rjv = 0.0) the flow con- 
straint is fulfilled but no stable hybrid stars are obtained. 



mimick the DBHF-EoS and generalize to finite tempera- 
tures (DD-F4). Fig. 6 shows the resulting phase diagram 
including the transition from nuclear to quark matter 
{vd = 1.030, rjv = 0.50) which exhibits almost a crossover 
transition with a negligibly small coexistence region and a 
tiny density jump. At temperatures beyond T 45 MeV 
our NM description is not reliable any more since contri- 
butions from mesons, hyperons and nuclear resonances are 
missing. This will be amended in future studies. 

5. Conclusions 

We have presented a new scheme for testing nuclear mat- 
ter equations of state at supernuclear densities using con- 
straints from neutron star and HIC phenomenology. Mod- 
ern constraints from the mass and mass-radius-relation 
measurements require stiff EsoS at high densities, whereas 
flow data from heavy-ion collisions seem to disfavor too 
stiff behavior of the EoS. As a compromise we have pre- 
sented a hybrid EoS with a phase transition to color su- 
perconducting quark matter which, due to a vector meson 
meanfield, is stiff enough at high densities to allow com- 
pact stars with a mass of 2 Mq. Such a hybrid EoS could 
be superior to a purely hadronic one as it allows a faster 
cooling of objects within the typical CS mass region. This 
way, young nearby X-ray dim objects such as Vela could be 
explained with masses not exceeding 1.5 Mq. The present 
hybrid EoS predicts hybrid stars that "masquerade" as 
neutron stars, suggesting only a tiny density jump at the 
phase transition. This characteristics is also present for 
the symmetric matter case and persists at higher temper- 



atures in the QCD phase diagram. It is suggested that the 
CBM experiment at FAIR might softly enter the quark 
matter domain without extraordinary hydrodynamical ef- 
fects from the deconfinement transition. 
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